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ABSTRACT 

It has long been thought that the pulsed X-ray properties of rotation-powered pulsars are stable 
on long time scales. However, long-term, systematic studies of individual sources have been lacking. 
Furthermore, dramatic X-ray variability has now been observed from two pulsars having inferred sub- 
critical dipole magnetic fields. Here we present an analysis of the long-term pulsed X-ray properties 
of the young, energetic pulsar PSR B1509— 58 using data from the Rossi X-ray Timing Explorer. We 
measured the 2-50 keV pulsed flux for 14.7 yr of X-ray observations and found that it is consistent 
with being constant on all relevant time scales, and place a 3cr upper limit on day-to-weck variability 
of <28%. In addition, we searched for magnetar-like X-ray bursts in all observations and found none, 
which we use to constrain the measurable burst rate to less than one per 750 ks of observations. We 
also searched for variability in the pulse profile and found that it is consistent with being stable on 
time scales of days to decades. This supports the hypothesis that X-ray properties of rotation-powered 
X-ray pulsars can be stable on decade-long time scales. In addition, we extend the existing timing 
solution by 7.1 yr to a total of 28.4 yr and report updated values of the braking index, n = 2.832±0.003 
and the second braking index, m — 17.6 ± 1.9. 

Subject headings: pulsars: individual (PSR B1509— 58), supernovae: individual (G320.4— 1.2) 



1. INTRODUCTION 

The X-ray emission properties of non-accreting neu- 
tron stars are today known to be extremely diverse. 
There are multiple classes of neutron stars identified - 
rotation-powered pulsars (RPPs), Anomalous X-ray Pul- 
sars (AXPs), Soft Gamma Repeaters (SGRs), thermally 
emitting isolated neutron stars, central compact objects 
- to name a few (see Kaspi 2010, for a review), with the 
basis for classification of a given neutron star into these 
different classes often being its X-ray emission properties. 

Anomalous X-ray Pulsars (AXPs) and Soft Gamma 
Repeaters (SGRs) are now well established as being ex- 
tremely variable in the X-ray range, showing a wide va- 
riety of often dramatic bursting behavior during active 
periods, in addition to lower-level flux and profile vari- 
ability during periods of relative quiescence (see Woods 
& Thompson 2006, Kaspi 2007, Mereghetti 2008, and 
Rea & Esposito 2011, for reviews). This variability is 
ascribed to activity both within the stellar interior as 
1 well as in the magnetosphere, both owing ultimately 
to the decay of these objects' extremely high magnetic 
fields {B > Bqed = 4.4 x 10^^ G), particularly in 
the context of the now wide l y accepted magnetar rn o del 
(jThompson fc DuncMlll996l : [Thompson et al.ll2002[ ). 

Interestingly, two recent results have lowered the esti- 
mated dipolar field strength necessary for magnetar-like 
activ ity to well below B ~ 10^** G (e.g. [Thompson et al.l 
[200l . The young RPP PSR J1846-0258 with B = 
4.9 X 10^'^ G experienced a magnetar-like outbu rst in 2006 
(|Gavriil et al.ll2008HKmmr fc Safi-Harbll2008D . while the 
SGR 0418-^5729, discovered via an outburst in 2009, has 
an upper Umit o n its dipole field of i? < 7.5 x 10^^ G 
(jRea et al.l 120101 ). Thus, the smallest possible B field 
needed to power a magnetar-like outburst is not cur- 

^ maggieQphysics.mcgill.ca 



rently known. These recent discoveries highlight the fact 
that the link between magnetars and RPPs is poorly un- 
derstood. In particular, it suggests that the X-ray emis- 
sion of some RPPs may be more variable than pre viously 
thought. Indeed, a study bv lPerna fc PonsI ((20ll argues 
that bursting behavior may be generically possible in all 
RPPs, only more common in those that are young and 
with high B. 

Systematic examinations of the X-ray variability prop- 
erties of RPPs in order to look for other examples of 
bursting behavior or fiux changes has been hampered by 
the difficulty in obtaining a well sampled data set with a 
single, well calibrated instrument over a long time span. 
Rotation-powered pulsars are generally too faint for de- 
tection with all-sky monitors. 

PSR B1509— 58 is a young, energetic pulsar 
that was discovered with the Einstein satellite 
(jSeward fc Harnden Jr.|[T98^ . It has period 150 ms, a 
high spin-down luminosity of ii^ = 1.7 x 10'^'' ergs~^, 
and a relatively large inferrecQ dipolar magnetic field 
of i? = 1.5 X 10^"^ G, larger than that for the mag- 
netar with the s mallest estimate of i? < 7.5 x 10^^ G 
(|Rea et al.l [20Tol ). PSR B1509-58's high spin-down 
rate and lack of glitches have allowed the measure- 
ment of higher-order spin parameters, in particular, 
the second and third frequency derivatives {i) and 
v). These provide a determination of the braking 
index, n = vv jiP' = 2.839 ± 0.003, and second braking 
index, m = v^ 'vjv^ = 18.3 ± 2.9 (|Kaspi et al.l 119941 : 
iLivingsto ne et al.ll2005b[ ). all owing tests of th e standard 
pulsar spin-down model fe.g. lBlandfordl[T994f) . 

PSR B1509— 58 presents an excellent opportunity for 
a systematic, long-term study of a RPP's X-ray emission 

^ Note that the dipolar magnetic field for pulsars is _B = 3.2 X 
10i'-*yP?G. 
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properties, since it has been observed regularly since 1996 
with the Proportional Counter Array (PCA) aboard the 
Rossi X-ray Timing Explorer (RXTE) as a timing cali- 
bration source. Here we present a search of 14.7 yr of 
RXTE data of PSR B1509-58 for pulsed flux variations, 
pulse profile variations, and magnetar-like bursts. We 
also extend the timing solution for the pulsar by 7.1 yr, 
increasing the total observing time to 28.4 yr. We use 
these new data to refine the measurements of n and m, 
and further quantify the variability of n. 

2. OBSERVATIONS 

PSR B1509-58 data were obta ined with th e Pro- 
portional Counter Array (PCA; iJahoda et al.l 119961 : 
iJahoda et all 12006^ on board RXTE. The PCA con- 
sists of five collimated xenon/methane multi-anode pro- 
portional counter units (PCUs) operating in the 2- 
60 keV range, with a total effective area of approxi- 
mately 6500 cm^ and a field of view of - 1° FWHM. 
Wc downloaded 290 archival RXTE observations of 
PSR B1509~58 from the HEASARC archiveQ. All data 
were collected in the GoodXenon mode, which records 
the arrival time (1 ^s resolution) and energy (256 chan- 
nel resolution) of each unrejected event. The observa- 
tions span 7.1 yr from MJD 50148-55221 (1996 March 6 
-2010 November 21). 

The data reduction and preparation is significantly dif- 
ferent for the timing and pulse profile analyses than for 
the pulsed flux analysis. Details for each are included in 
the following sections. 

3. TIMING ANALYSIS 

Our timing analysis of PSR B1509— 58 follows the com- 
mon phase-coherent approach, in which we account for 
each rotation of the pulsar. Pulse times of arrival are 
fitted with a Taylor expansion of pulse phase (see, for ex- 
ample, iManchester fc Tavlorili977f ). allowing the extrac- 
tion of precise spin parameters (spin-frequency and 
frequency derivatives t', i', and v). 

In order to produce pulse times-of-arrival (TOAs), we 
extracted photons from all three xenon layers of each 
PCU in the 2 - 30keV energy range (or 4-72 channel 
range). Data from the individual PCUs were merged 
and binned at 1/1024 s resolution. The time series 
were reduced to barycentric dynamical time (TDB) at 
the solar system barycenter using the known position 
from radio interferom etry (J2000 RA = I S'' 13'"55';62, 
Dec = -59° 08' 9.00"; IGaensler et aII[T999h and the JPL 
DE200 solar system ephemeris. 

When more than one observation occurred in a 24- 
hr period, we combined the data from multiple observa- 
tions into a single time series in order to produce better 
TOAs and having smaller uncertainties. Each time series 
was f olded with 64 phas e bins u sing an initial ephemeris 
from [Livingstone et al.l (|2005bf ). Resulting pulse pro- 
files were cross-correlated in the Fourier domain with 
a high signal-to-noise-ratio template. We implemented 
a Fourier domain filter by using six harmonics in the 
cross-correlation procedure. Cross-correlation produces 
an average TOA for each observation with a typical un- 
certainty of '^0.5 ms. TOAs were fitted phase coherently 
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Figure 1. RXTE timing residuals for PSR B1509— 58 spanning 
1996 March - 2010 November. The top panel has v, v, and u fitted; 
the middle panel has a quartic polynomial (i.e. v fitted), and the 
bottom panel has 12 derivatives fitted. Spin parameters from a 
fully phase-coherent fit are given in Table m 

with the pulse timing software package TEMPC0; timing 
residuals and refined spin parameters are produced as 
output. Spin parameters for all 14.7 yr of X-ray data are 
given in Table [TJ 

The top panel of Figure [1] shows timing residuals with 
pulse frequency, frequency derivative, and second 
frequency derivative, v, fitted. No sudden jumps were 
found during the phase connection process, nor were 
any cusps identified in the timing residuals, indicating 
that no glitches occurred during the time spanned by 
the RXTE observations. The middle panel shows resid- 
uals with ii also fitted. PSR B1509— 58 is one of two 
pulsars with deterministic measurements o f v and v 
(jKaspi et al.l 119941: [Livingstone et al.l l2005bD . however, 
after the removal of all deterministic parameters, sig- 
nificant structure remains in the timing residuals, as is 
visible in the middle panel. This can be attributed to 
timing noise, a low-frequency or "red" noise process that 
is a common but poorly understood phenomenon in pul- 
sars (and previously disc ussed in detail for this pulsar in 
iLivingstone et al.|[2005b[) . TEMPO has the capability to fit 
up to 12 frequency derivatives; the results of this fit are 
shown in the bottom panel of Figure [TJ The presence 
of timing noise is known to contaminate measurements 
of higher order spin parameters such as n and m and 
complicates the proper estimation of uncertainties. One 
method of measuring determinsitic parameters and esti- 
mating their true uncertainites in the presence of timing 
noise is discussed in the following section. 

3.1. Partially phase- coherent timing analysis 

One method of finding reliable values for n and m 
in the presence of timing noise is to use a piece-wise 
approach, deemed "partially coherent" timing analysis. 
The data are broken into subsets, each of which is fitted 
phase-coherently with v, v, and v. For each subset, we 
require that no structure be visible in the resulting timing 



^ http:/ /heasarc. gsfc.nasa.gov/docs/archivc. html 
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residuals. The amount of data included in each subset 
is determined only by this requirement, so is allowed to 
vary. On average, the time span in a subset is '^SOOdays. 
The date ranges spanned by these subsets and the result- 
ing ephemerides are given in Table [2] (and are used for 
our pulsed flux analysis, see Section|4|). In order to better 
sample the available data, we created additional subsets 
that overlapped the original subsets by ~ 1/2. Overlap- 
ping subsets are useful since short coherent timing so- 
lutions produce fitted parameters dominated by the end 
points (e.g. Livingstone et al. 2005b). Figure [5] shows 
all measurements of v both from the current analysis, 
as well as previous measurements determined using the 
same meth od but from both radio and X-ray timing data 
reported in [Livingstone et all ()2005b[ ). 

We performed a weighted least squares fit to mea- 
sure i) and a bootstrap error analysis to better estimate 
the uncertainties, since the data are known to be con- 
taminated with red noise leading to formal uncertaintie s 
that underestimate the true uncertainties ()EfronllT979[ ). 
This analysis results in i>' = -1.21(13) x 10"^^ s"", in 
agreement with the previously reported value of i) = 
-1.28(21) X 10-31 s"^, based on 21.3yr rather than 
28.4 yr of data, but calculated in the same way. 

A value of can be used to find a measurement of the 
"second" braking index, m, via m = v^ 'v jv^. This pa- 
rameter, a higher order analog to the braking index, n, 
provides a unique test of p ulsar spin-down models (e.g. 
iBlandford fc Romanil ITosl ) . The above value of V cor- 
responds to a value of m = 17.6 ± 1.9 . smaller than the 
value fr om 21.3 yr of data reported in [Livingstone et al.l 
(|2005bf) of TO = 18.3 ± 2.9, but not significantly so. 

We used the same timing solutions to calculate n for 
each data subset as shown in Figure [2] The weighted 
average value, with the uncertainty calculated from a 
bootstrap analysis is n = 2.832 ± 0.003, smaller than the 
previously reported value of 2.839 ± 0.003 at the 1.6cr 
level. In order to look for a linear change in n over the 
entire 28.4 yr of observations, we performed a weighted 
linear least squares fit to the data. This results in a slope 
of An/Af = (— 4±2) x 10"^, which is consistent with no 
secular change in n at the 2ct level. While the long-term 
n is relatively stable, short-term variability is visible as 
scatter in the individual measurements, which can most 
likely be attributed to timing noise, since no discrete 
jumps in i.e. glitches, were found. This amounts to a 
maximum variation from the average value of 4.8 ±0.3%, 
and a root-mean-square variation from the mean of 3.5%. 
Because of the large scatter in the measurement of v near 
MJD 55000, which is significantly larger than the mean, 
we further verified that no glitch occurred in this time 
period. 

4. PULSED FLUX ANALYSIS AND RESULTS 

In order to determine the pulsed fiux of PSR B1509— 58 
at each epoch, for each observation, we created a phase 
resolved spectrum with 8 phase bins, for each PCU. 
Each observation was folded with the software package 
fasebir{3 using the appropriate short coherent ephemeris 
described in Section 13.11 and listed in Table [2] We 
used these short ephemerides instead of the global so- 
lution in part because fasebin can take a maximum 
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Figure 2. Braking indices from short phase- coherent timing so- 
lutions. Data prior to MJD 53000 were previously reported in 
jLivingstonc ct al. 2005b). The weighted average is n = 2.832 it 
0.003. At the 2(T level, there is no linear change over 28.4 yr. Signifi- 
cant scatter is apparent in the data; the root-mean square variation 
from the mean n is at the 3.5%-level. 

of two frequency derivatives for an input ephemeris. 
For PSR B1509— 58, using the global two-derivative 
ephemeris produces very poor pulse profiles because 
of the effect of the deterministic ii and timing noise, 
whereas the profiles produced from short ephemerides 
are consistently good. A phase-averaged spectrum was 
also created in order to build a response matrix for each 
PCU. For each observation, we combined phase-resolved 
spectra and response matrices for all PCUs that were on 
for > 95% of the observation. 

From each combined phase-resolved spectrum, we cre- 
ated a pulse profile with 8 phase bins in the 2-50 keV 
energy range. The minimum bin of each profile was used 
to set the background level, and was subtracted from all 
bins of the phase-resolved profile. 

After background subtraction, we combined the spec- 
tra for all phase bins into a single spectrum and fit- 
ted the result using XSPECQ. For each observation, we 
used a simple power-law model with a fixed spectral 
index of F = 1 19 an d fixed = 8.6 x lO^^ cm^^ 
(jCusumano et al.]|200lD . fitting only for the flux. The 
resulting fit values showed that this model fit all of 
our data well. 

For short exposure times, there is an observed anti- 
correlation between the measured pulsed fiux and the 
source exposure time. This is a statistical effect present 
for all sources and disappears when sufficient counts are 
detected, i.e. for longer exposure times. While one might 
expect that the uncertainties on pulsed flux measure- 
ments would simply increase for short exposures, in ad- 
dition, the measured pulsed fiux is biased upwards. This 
arises from Poissonian variations, which lead to the mini- 
mum bin being biased downward when there are few total 
counts, thus to an under-estimate of the background and 
an over-estimate of the pulsed fiux. Because the number 
of active PCUs changes from observation to observation, 
for convenience, we define an "effective" observation time 
of Tcff = A^pcusTobs, where N^^ms is the number of PCUs 
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Figure 3. Absorbed pulsed flux measurements for PSR B1509— 58 
determined for observations with T^ff > 6ks. Vertical lines mark 
dates of epoch changes in RXTE. Filled circles (colored red in 
the online version) show the average pulsed flux over year-long 
periods. The average pulsed flux in the 2-50 keV energy range 
is 1.181(5) X 10~^^ ergcm"^ s~^. A weighted least-squares flt to 
the data gives a slope of (— 7.4 ± 3.8) X 10~^^ ergcm"^ day~^ , 
indicating that the pulsed flux is consistent with having no linear 
trend at the l.9cr level. 

on during an observation and Tobs is the on-source time. 
Wc compared the measured flux values to Toff and found 
a small bias to be present. This bias was eliminated 
by excluding all observations having T^ff < 6 ks, which 
amounts to 34% of the data. 

Figure [3] shows all measurements of the pulsed com- 
ponent of the flux of PSR B1509— 58 for observations 
with Tcff > 6ks. Vertical lines delimit each of the RXTE 
epochs (1, 2, 3, 3a, 4) and the losses of the propane layers 
in PCUO (also referred to as epoch 5) and PCUlO]- We 
verified that there are no significant correlations between 
the pulsed flux and any of the gain changes or propane 
layer loss events. We also measured the pulsed count 
rate for each observation and note that correlations ex- 
ist between the pulsed count rate and gain changes and 
propane layer losses. In addition, the pulsed count rate 
decreases gradually by ~10% over the lifetime of RXTE, 
associated with the long-term degredation of the PCA. 
For these reasons, we report only the pulsed flux val- 
ues, despite somewhat larger uncertainties. We note that 
PSR B1509— 58 is listed as a calibration source for RXTE 
but has only been used as a timing calibrator and not a 
flux calibration source^- 

The weighted average 2-50 keV absorbed pulsed flux 
for all 14.7yr is 1.181(5) x 10"^" ergcm"^ s~^, with a 
value of 52.4 for 150 degrees of freedom. A weighted 
linear least-squares fit to the pulsed flux values / gives 
Af/At= (-7.4 ±3.8) X 10-16 erg cm-^s-^day-i, that 
is, consistent with having no linear trend at the 1.9a 
level. The maximum variation of the pulsed flux of a sin- 
gle observation from the average corresponds to a vari- 
ability of (10±6)%. At the 3(7 level, we can rule out 
changes on day-to-week time scales of greater than 28% 
from the average pulsed flux value. We also studied the 

^ http: / /heasarc. gsfc.nasa.gov/docs/xte/pca _history.html 

|http: / /heasarc. gsfc.nasa.gov/docs/xte/e-c_table. html 

* See, for example, http://heasarc.nasa.gov/docs/xte/abc/time.html 
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Figure 4. Summed 64-bin pulse profile for PSR B1509— 58 for all 
14 Cycles of RXTE data for the 2-50 keV energy range. We found 
no significant difference between individual pulse profiles and the 
average template, as discussed in Section [5] 

flux stability on monthly and yearly time scales by creat- 
ing weighted averages of individual flux measurements. 
Average flux values for each RXTE cycle are shown on 
Figure [3] as filled circles (colored red in the online ver- 
sion) On monthly time scales, at the Za level, we rule 
out variation from the average larger than 20%, while on 
yearly time scales, the 3cr upper limit on variability is 
21%. 

5. PULSE PROFILE ANALYSIS 

In order to search for variability in the X-ray pulse 
profile of PSR B1509-58, we folded each 2-30 keV time 
series (as described for the timing analysis in Section [3]) 
with the appropriate ephemeris given in Table [5J Pre- 
vious studies have examined the pulse profile as a func- 
tion of energy in detail and found that the pulse shape 
does not ch ange appreci ably in t he RXTE energy range 
(jRots et al.1.1998: Marsden et al.l il997') , hence we do not 
repeat this analysis here and concentrate instead on the 
possibility of the profile changing with time. 

First, we analyzed the profiles for individual observa- 
tions in order to search for variability on time scales of 
days to months (i.e. the range of time spans between 
individual observations). We created a template profile 
by summing all aligned pulse profiles, shown in Figured 
We subtracted a DC offset from both the individual pro- 
files and template, then normalized the template to each 
individual profile, subtracted them and found the resid- 
uals between the two profiles. We calculated the reduced 
of the resulting residuals and in each case, we see no 
evidence for pulse profile variability. To quantify the lack 
of profile changes, we calculated the root-mean-square of 
the profile residuals as a percentage of the number of 
counts in each DC-ofFset subtracted profile and found a 
3cr upper limit on variations of 2.1%. 

Next, we created a single profile for each RXTE ob- 
serving cycle (spanning ~l-yr each). In order to check 
whether there were any small-scale changes in the pulse 
profile, we compared each profile against the previously 
described template profile. Repeating the analysis de- 
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scribed above, again we found no evidence for variabil- 
ity. For the yearly profiles, we found a upper limit on 
profile changes of less than 0.1%. 

We also compared the profile from each of the 14 
RXTE cycles with that of the previous and subsequent 
cycle profiles. We normalized each profile and then cre- 
ated residuals between pairs of profiles. We then per- 
formed a test on each set of residuals. In every case, 
< 1, i.e. we detected no evidence of long term 
changes in the pulse profile. In order to have maximum 
sensitivity to changes on the longest possible time scales, 
we also repeated the analysis by comparing the Cycle 1 
profile with the Cycle 14 profile. Again, we found that 
X\ < 1. 

6. BURST SEARCH 

To search for magnetar-like X-ray bursts from 
PSR B1509— 58, for each observation we created a time 
series for each active PCU from GoodXenon data, se- 
lecting events in the 2 - 20keV range and from all 
three detection layers (the same energy range as selected 
for sim ilar searches for X-ray bursts from AXPs, e.g. 
iGavrii l et al. 2004). This was repeated several times 
with different time resolutions, specifically for 15.625 ms, 
31.25 ms, 62.5 ms and 125 ms, to provide sensitivity to 
bursts on a hierarchy of time scales. The procedure use d 
here was based on that described in'Gavr iil et al.l ()2004[ ). 

By this algorithm, bursts are identified by comparing 
the count rate in the i*'' bin to the average count rate. 
Because the background rate of the PCA typically varies 
over a single observation, we calculated a local mean 
around the i*'' bin. We then compared the count rate 
in the i*'* bin to the local average. If the count rate in 
a single bin was larger than the local average, the prob- 
ability of such a count rate occuring by chance was cal- 
culated. Then the probability of the count rate in the 
corresponding bin in the other active PCUs was also cal- 
culated (whether or not the count rate in that bin was 
greater than the local average) . If a PCU was off during 
the bin of interest, its probability was set to 1. We then 
found the total probability that a burst was observed. 
If the total probability of an event was P^^tot < 0.01/iV 
(where N is the total number of time bins searched) , it 
was flagged as a burst. 

We searched a total of 751 ks of RXTE data using the 
above described process and found no X-ray bursts from 
the direction of PSR B1509-58. 

Unlike the AXPs, PSR B1509-58 has a bright PWN 
which contributes a large unpulsed background to the 
1° field-of-view of the PCA. Because of this larger back- 
ground, the sensitivity to AXP-like X-ray bursts could 
be significantly smaller than for the AXPs. Thus, in or- 
der to determine the fiuence of bursts that would have 
been detectable in these data, we added simulated bursts 
to real RXTE time series of PSR B1509-58. Each simu- 
lated burst was confined to a single time bin and spread 
between the operational PCUs. We added many bursts 
to each time series, each separated by sufficient time that 
only one burst would ever be included in a calculation 
of the running mean. Many bursts were added to each 
time series in order to include the effect of the variation 
of the background count rate, largely owing to the vari- 
able particle background during the RXTE orbit. We 
repeated this analysis for each time resolution searched. 



Unsurprisingly, we found that the minimum burst fiu- 
ence required to detect 100% of the simulated bursts is 
a function of both the number of operational PCUs and 
width of the time bins. 

In terms of detected burst count rates, the smallest 
detectable bursts are similar to the smallest bursts de- 
tected from the AXP IE 2 259-^586 in its large 2002 out- 
burst ijCavriil et al.l [200^ . For example, for 31.25 ms 
time resolution, 100% of simulated bursts are detected 
if the burst consisted of 8 or more counts per PCU 
or more than 25 counts/burst total (i.e. at least 25 
counts per 31.35ms time bin). In order to quantify 
this in terms of burst fiuences, we assumed a burst 
spectrum of a simple power law with spectral index 
of r ~ 1.35, as measured for the 80 bursts found 
from the AXP IE 2259-f586 (|Gavriil et al.l[200l . For 
62.5 ms time resolutions, we found that burst fiuences of 
> 1.5 X 10~^^ergcm~^ (2-20 keV) would have been de- 
tected, an order of magnitude less than the smallest burst 
fiuenc e detected from PSR J1846-0258 in its 2006 out- 
burst (jCavrnl et al.ll2008h . Detectable burst rates and 
fiuences for all searched time resolutions are given in 
Table [Sj and are valid in the absence of instrumental 
flares. Thus, we can confidently exclude the possibil- 
ity that AXP-like X-ray bursts with fiuences larger than 
the limits given in Table [3] were present in the RXTE 
PSR B1509-58 data. 

7. DISCUSSION 

We searched for several types of X-ray variability from 
the young, high-i? field pulsar PSR B1509— 58 and found 
that the pulsar is stable in each way studied. This 
establishes a baseline for the quiescent pulsed flux of 
PSR B1509— 58 in case of a future magnetar-like out- 
burst. We establish a 3f7 upper limit on pulsed flux 
variations of less than 28% from the average, ruling out 
pulsed flux changes similar to those seen in the magne- 
tars and PSR J1846-0258. For example, during the 2006 
outburst of PSR J1846-0258, the 2-60 keV pulsed flux 
as measured with R XTE increased by a factor of ^ 3 
(|Gavrnl et al.l [200g). 

The burst search algorithm that has successfully 
discovered X-ray bursts from several AXPs and 
PSR J1846— 0258 did not reveal any bursts from the di- 
rection of PSR B1509-58. We searched 751 ks of RXTE 
observations for X-ray bursts, implying that the burst 
rate must be less than one burst per 751 ks. This is less 
than the burst rate for PSR J1846-0258, where 5 X- 
ray bursts hav e been observed, implyi ng a burst rate of 
one per 407 ks (jLivingstone et al.ll2011h . Bona fide AXPs 
have a wide range of observed burst rates, from < 1 burst 
per -700 ks in IRXS J170849. 0-400910, up to 1 burst 
per -30 ks for IE 2259+586. 

Finally, a study of the X-ray pulse proflle of 
PSR B1509— 58 shows no significant variability on time 
scales of days to years. Variations in individual pulse 
profiles are constrained to be less than 2.1% from the 
average profile at the 3(7 level. Together, these results 
confirm the hypothesis that RPPs can have stable X-ray 
properties on both short- and long-term time scales. 

7.1. Radiative variability from rotation-powered pulsars 

^ Previous studies have shown a correlation between 
E and Lx in the soft X-ray band for magneto- 
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spheri c emission from RPP s and their nebular emis- 
sion (iSeward fc Wand 119881: iBecker fc Trumped 119971: 
iKargaltsev k Pavlovll2008i r This can be interpreted as a 
constant value of the conversion of spin-down luminosity 
into X-ray luminosity, rjx = Lx /E, for all RPPs. Estab- 
lishing a similar correlation between pulsed X-ray lumi- 
nosity and E would be interesting, and implicitly requires 
Lx to be stable. Prior to this study of PSR B1509-58, 
the long-term stability of Lx had not been well estab- 
lished for any pulsar, largely owing to systematic uncer- 
tainties when comparing results from different telescopes. 
Here, we show that Lx can be steady on day-to-year-long 
time scales for RPPs. This confirms the initial hypothe- 
sis, and stands in stark contrast to the magnetars, where 
Lx varies dramatically. 

Recent studies have shown that radio emission vari- 
ability is m ore common in pulsars than previousl y 
thought (e.g. iLvne eralll2010t IWeltevrede et all 1201 If ). 
Specifically, the radio emission varia.bility seen in the 
"intermittent" pulsars (iKramer et all 120061: iLvne et al.l 
[2010) appears to be associated with large-scale magne- 
tospheric variations. If such a mechanism were active in 
PSR B1509— 58, variations in Lx would be expected, and 
should be visible in our data, along with variations in the 
X-ray pulse profile. Because high-energy pulsar emission 
is more closely associated with the outer magnetosphere, 
naively, one might expect more dramatic variations in the 
X-ray regime than in the energetically unimportant radio 
regime. So far this cannot be confirmed, because none 
of the "intermittent" pulsars have detectable magneto- 
spheric X-ray pulsations, and no variations in the X-ray 
pulse profile of PSR B1509— 58 were detected here. 

Interestingly, two high B-field RPPs, PSR J1846-0258 
and J1119— 6127, have shown unusual radiative behav- 
ior near the time of a rotational glitch, in the X- 
ray and radio regimes, respect i vely fe uiper & Hermsen 
[200l [Livingstone eTall [2010l : IWeltevrede et al.l 1201 lO . 
Given this apparent propensity for pulsars with higher- 
than-average i?-fields to display radiative variability and 
PSR B1509— 58's relatively large magnetic field of 1.5 x 
10^'^ G, future observations of both the radio and X-ray 
emission properties of PSR B1509— 58 may show unusual 
behavior, particular in the case of a glitch. 

7.2. Relationship between B and magnetar outbursts 

The exact relationship between the estimated dipole 
field strength and the ability and/or frequency of a neu- 
tron star to exhibit magnetar activity remains to be es- 
tablished. However, there is now a wide range of dipole 
B estimates in neutron stars exhibiting magnetar out- 
b ursts. 

iPerna fc Pons! ()2011f ) suggest that all neutron stars 
may exhibit the type of X-ray outbursts typically associ- 
ated with magnetars, but that the frequency of events is 
highly dependent on both the age and the magnetic field 
strength of the source. For very young neutron stars 
(<2kyr), outbursts may be observed roughly yearly for 
B ~ 8 X lO^^^G, but only every 50-100 yr for B ~ 2 x 
lO^'' G. Thus for neutron stars with B ~ 10^^ G, such as 
PSR B1509-58 and PSR J1846-0258, outbursts should 
be less frequent than once per century. This implies that 
the observation of an outburst from PSR J1846— 0258 
was relatively unlikely given that it has only been ob- 
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Figure 5. Measurements of i> from short phase-coherent timing 
solutions. Data prior to MJD 53000 were previously reported in 
ILivingstone et ahl l l2005bl ). A weighted least-squares fit to the data 
gives a slope of i> = 1.21(13) X 10~^^ s~*. The effect of timing 
noise is visible as scatter in the data. 

served for '^11 yr, and that the lack of magnetic activity 
from PSR B1509— 58 is not surprising. Ongoing obser- 
vations of many high-B radio pulsars should increase the 
chances of further observations of magnetar-like bursts 
from RPPs. 

7.3. Braking index, timing noise, and glitches 

We presented an extention of the timing solution o f 
PSR B1509-58 from 21.3 yr (jLivingstone et al.l [20050 ) 
to 28.4yr. We found a braking index of n = 2.832(3), 
which is 1.6(7 from the value based on 21. Syr of data. 

The scatter in short measurements of i), and therefore 
in 71, is visible in both Figures [5] and [2l and is typical of 
timing noise in young pulsars. We observed that the root- 
mean-square variation in n is 3.5%, larger tha n the varia- 
tions observed in the Crab pulsar of ~ 0.5% (jLvne et al.l 
119931 ). We note that there is evidence for timing noise, 
but no eviden ce for a glitc h in c ontrast to the sugges- 
tion made by iHobbs et al.l (|2010[ ) that timing noise in 
pulsars with characteristic ages less than Tc < 10^ yr is 
entirely because of unmodeled glitch recovery. However, 
it remains possible that PSR B1509— 58 experiences very 
small glitches that cannot be distinguished from timing 
noise given the current data. 

Because no glitches were found in the latest 7.1 yr of 
timing data, this extends the time without a detected 
glitch for PSR B1509-58 to 28.4 yr. While the youngest 
pulsars (tc < 2kyr) tend to have smaller and fewer 
glitches than somewhat older pulsars (tc ^ 5 — lOkyr), 
PSR B1509— 58 is extreme in its lack of glitches. For ex- 
ample, the similarly aged Crab pulsar has had 24 glitches 
in 42 yr of observations, PSR B0540— 69 has had one 
glitch in 7.6 yr, and PSR J1119-6127 has had 3 glitches 
in 12 vr (lEspinoza et al.ll20ll[: ILivingstone et al.|[2005al: 
IWeltevrede et al.ll2011h . 

A higher order analog to the braking index (n = 
vvjiJ^^, is the "second braking index", which can be pre- 
dicted from a measurement of n to be mg = n(2n— 1) and 
measured with a value of v to be m — iiv^ jir" . Timing 
noise and glitches prevent a measurement of n for most 
pulsars. A measurement of m is even more rare - signif- 
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icant measurements have only been rnade for the Crab 
pulsa r and PSR B1509-58 (iLvne et al.ll993t iKaspi et alj 
Il994t iLivingstone et all l2005b[ ). If the standard spin- 
down model {i^ = —Kv'^) is an accurate description of 
the spin down of the pulsar, then m — utiq. By contrast, 
n < 3 and m ^ niQ implies that K = K{t) and a mea- 
surement of m different fro m mp provides constraints o n 
the functional form of K{t) (|Blandford fc Romanilll988D . 

Using the full span of available data of 28.4 yr, we found 
ly' = -1.21(13) X 10-31 s"'', which is within la from 
the previously reported value based on 21.3 yr of data. 
For the current measurements of n and v , we find that 
mo = 13.21 ± 0.03 and m = 17.6 ± 1.9, which differ at 
the 2.3cr level (see Table[T]for all rotational parameters). 
If uncertainties continue to decrease at the same rate as 
over the last 7.1 yr since the previous measurement, and 
the measured value of m remains stable, with another 
^Qyv of timing observations, the uncertainties should be 
small enough that m would differ from toq at the 3ct level. 
If this proves to be the case and the two values are in- 
consistent, this would present chall enges to the standard 
model of pulsar spin-down (e.g., IBlandford fc Romanil 
[1988; ,Blandford.l99 4l . 

8. CONCLUSIONS 

Analyzing 14.7yr of RXTE data, we found no evidence 
of variability in the X-ray emission of PSR B1509— 58. 
We have placed limits on possible changes in the 
pulsed flux, profile variability, and burst rate from 
PSR B1509— 58. This indicates that it is indeed possible, 
as typically assumed, for the X-ray emission properties 
of RPPs to be stable on day to decade long time scales. 

Several RPPs have dipole B estimates close to, or 
larger than, magnetar-strength fields. These RPPs are 
typically referred to as the high B-field pulsars. The 
link between these and magnetars is not well under- 
stood and is und er investigat ion (e.g. Olausen et al.ll2010l : 
IZhu et al.l I20T1I |2009; Ng fc Kaspil I2010l ). These pul- 
sars, along with PSR J184 6— 0258, are possibly relate d 
to the transient magnetars (IKaspj fc McLaughlin|[2005l) . 
such as XTE 1810—197 which experience long periods of 
quies cence, followed by magnetar activity ([Ibrahim et al.l 

If the internal B field of a neutron star can be much 
larger than the observed external field responsi ble for the 
sgin-down (as suggested for SGR 0418+5729: iRea et al.l 
[2OIO1) then it becomes difficult to predict which pulsars 
may experience a future magnetar-like outburst. Indeed, 
long-term monitoring observations of PSR J1846— 0258 
showed that a sudden and dramatic (though temporary) 
change from a pulsar to a magnetar is possible. Thus 
X-ray monitoring of pulsars such as that performed with 
RXTE of PSR B1509-58 provides an opportunity to dis- 
cover future outbursts from less-than-obvious magnetar 
candidates. 
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Table 1 

X-ray Spin Parameters for PSR B1509-58 



Parameters for phase-coherent analysis 



Dates (Modified Julian Day) 
Epoch (Modified Julian Day) 
Ntoa 

V (Hz) 

V (10-" s-2) 
v (10-21 s-») 
v (10-31 s-4) 



50148.096 - 55521.082 
52834.589000 
188 

6.611515243850(3) 
-6.694371307(6) 
1.9185594(5) 
-0.9139(2) 



Parameters for partially coherent analysis 



Average Braking Index, n 
mo = n{2n — 1) 



V (10- 



Second braking index, m 



2.832(3) 
13.21(3) 
-1.21(13) 
17.6(1.9) 



Note. — Spin parameters for 14.7 yr of X-ray timing data for PSR B1509— 58. 



Table 2 

Partially coherent ephemerides 



Epoch 


V 


V 


V 


Date range 


(MJD) 


(Hz) 


(10-11 Hz-2) 


(10-21 Hz-3) 


(MJD) 


50448 


6.6253600487(2) 


-6.734167(2) 


1.939(3) 


50148-50749 


51083 


6.6216683438(2) 


-6.7234900(8) 


1.932(2) 


50750-51449 


51853 


6.6171996051(2) 


-6.710680(1) 


1.917(3) 


51470-52167 


52616 


6.6127798918(2) 


-6.698028(1) 


1.912(3) 


52265-52926 


53091 


6.6100326310(1) 


-6.6901006(7) 


1.950(2) 


52893-53289 


53647 


6.60682103661(5) 


-6.6809030(3) 


1.9032(6) 


53311-53982 


54345 


6.60279542362(6) 


-6.6694306(4) 


1.9165(7) 


54007-54681 


54892 


6.5996455256(2) 


-6.660372(2) 


1.994(6) 


54707-55076 


55241 


6.5976380501(1) 


-6.6546543(5) 


1.894(1) 


54960-55522 



Note. — Uncertainties are formal TEMPO uncertainties. These spin parameters are used for the pulsed flux analysis; see Section |4] 



Table 3 

Minimum burst fluonccs to detect bursts 



Time resolution 


Total counts 


Counts/PCU 


Burst flucncc 


(ms) 






(10-12 erg cm-2) 


15.625 


20 


7 


2.6 


31.25 


25 


8 


5.9 


62.25 


30 


10 


15 


125.0 


40 


13 


39 



Note. — All fluences are for the 2— 20keV energy range. Fluences 
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